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ON ■ The production of antiprotons in proton-nucleus and nucleus-nucleus re- 

actions is calculated within the relativistic BUU approach employing proper 
selfenergies for the baryons and antiprotons and treating the p annihilation 
"q ' nonperturbatively. The differential cross section for the antiprotons is found 

to be very sensitive to the p selfenergy adopted. A detailed comparison with 
the available experimental data for p-nucleus and nucleus-nucleus reactions 
shows that the antiproton feels a moderately attractive mean-field at normal 
^ ■ nuclear matter density pq which is in line with a dispersive potential extracted 

^ . from the free annihilation cross section. 



The production of particles at energies below the free nucleon-nucleon thresh- 
old ('subthreshold production') constitutes one of the most promising sources of in- 
formation about the properties of nuclear matter at high densities since the particles 
are produced predominantly during the compressed stage at high density |], §, [3|. 
Antiproton production at energies of a few GeV/u is the most extreme subthreshold 
production process and has been observed in proton-nucleus collisions already more 
than 20 years ago @, |, §. Experiments at the JINR and at the BEVALAC 
H H have provided, furthermore, first measurements of subthreshold antiproton 
production in nucleus-nucleus collisions. Since then the problem was taken up again 
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at KEK JTDJ and GSI ]TTj] with new detector setups. Various descriptions for these 
data have been proposed. Based on thermal models it has been suggested that the 
antiproton yield contains large contributions from AiV — > p + X, A A — > p + X and 
PP ~ * production mechanisms [12|, [13|, |I3J| . Other models have attempted to 
explain these data in terms of multiparticle interactions Jl5] . 

In a first chance nucleon-nucleon collision model (assuming high momentum 



tails consistent with data on backward proton scattering) Shor et al. |L6] succeeded 
in reproducing the antiproton yield for the proton-nucleus case, however, underesti- 
mated the yield by more than 3 orders of magnitude for nucleus-nucleus collisions. 
This problem was partly resolved by Batko et al. [|l7j who performed the first 
nonequilibrium p-production study on the basis of the VUU transport equation. 
Within this approach it became possible to describe simultaneously the p-data from 
p + A and A + A reactions; however, the yield was still underestimated when in- 
cluding the strong p annihilation. Nevertheless, it became clear, that in A + A 
reactions the dominant production channel proceeds via an intermediate nucleon 
resonance which allows to store a sizeable amount of energy that can be used in a 
subsequent collision for the production of a pp pair. Later on, these results were 
also confirmed by Huang et al. |18|| within the QMD model; the authors achieved a 
reproduction of nucleus-nucleus data only when neglecting the annihilation channel. 

These results have led to the suggestion that the quasi-particle properties of 
the nucleons might be important for the p production process which become more 
significant with increasing nuclear density. Schaffher et al. [HJ found in a thermal 
relativistic model, assuming kinetic and chemical equilibrium, that the p-abundancy 
might be dramatically enhanced when assuming the antiproton selfenergy to be 
given by charge conjugation of the nucleon selfenergy. This leads to strong attrac- 
tive vector selfenergies of the antiproton. However, the above concept lacks unitarity 
between the real and the imaginary part of the p-selfenergy and thus remains ques- 
tionable. Besides this, even in the a — uj model the Fock terms lead to a suppression 
of the attractive p-field ||20|| , such that the production threshold is shifted up in 
energy again as compared to the simple model involving charge conjugation. Fur- 
thermore, the assumption of thermal and especially chemical equilibrium most likely 



is not fullfilled e.g. in Si + Si collisions around 2 GeV/u [21 



First fully relativistic transport calculations for antiproton production includ- 
ing p annihilation as well as the change of the quasi-particle properties in the medium 
have been reported in There it was found that according to the reduced nu- 
cleon mass in the medium the threshold for p-production is shifted to lower energy 
and the antiproton cross section prior to annihilation becomes enhanced for Si + 
Si at 2.1 GeV/u by roughly 2 orders of magnitude as compared to a relativistic 
cascade calculation where no in-medium effects are incorporated. The enhancement 
of the antiproton yield then was dramatically reduced again when including the 
strong annihilation channel which, however, lead - accidentally, as noted in [22| - to 
a reasonable reproduction of the p data for Si + Si at 1.65 and 2.1 GeV/u. A final 
answer on the antiproton problem especially with respect to the p selfenergy in the 
medium could not be given in since the explicit momentum dependence of the 
nucleons scalar and vector selfenergies had not been accounted for. 
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In this letter we therefore analyze the production of antiprotons within the 
framework of the relativistic transport theory (RBUU) that calculates p production 



perturbatively as in |T7|, |22| from the channels NN — > NNpp, AN — > NNpp as 
well as AA — > NNpp. The quasi-particle properties, i.e. the nucleon selfenergies 
U s (p, p),U fM (p, p), are taken from |24|, |25|] and are in line with the results of 



Dirac-Brueckner calculations [2~B| whereas the antiproton selfenergies are described 
on the basis of the a — u model [^7] with free coupling parameters and g%. 
Antiprotons are propagated explicitly in the respective time dependent potentials 
and their annihilation is calculated nonperturbatively by means of individual rate 
equations. A comparison with the most recent data from KEK and GSI on p- 
production will allow to approximately determine the free parameters gf and g% and 
provide first information on the antiproton potential in the medium. We finally 
compare the numerical results for the Schrodinger equivalent antiproton potential 
with a dispersive potential as evaluated from the annihilation rate and the charge 
conjugate potential as expected from simple relativistic mean-field theory. 

Since the covariant BUU approach has been extensively discussed in the re- 
views we only recall the basic equations and the corresponding quasi-particle 
properties that are required for a proper understanding of the results reported in this 
study. The relativistic BUU (RBUU) equation with momentum-dependent mean- 
fields or selfenergies is given by 

{[n,-n^fn + M*(^ (i) 

where f(x,p) is the Lorentz covariant phase-space distribution function, I co u is a 



collision term given in Ref. [pq| , and U s and are the scalar- and the vector- 

^tive mass M* and the kin 

terms of the fields by 



mean-fields. The effective mass M* and the kinetic momentum II M are defined in 



n M (x,p) = Pp-U^x^p) (2) 
M*(x,p) = M + U s (x,p), (3) 

while the quasi-particle mass-shell constraint is obtained as 

V(x,p)f(x,p) = (4) 

with the pseudo potential 

V(x,p) = l -{Ii\x,p) - M* 2 (x,p)}. (5) 

The above equation implies that the phase-space distribution function f(x,p) is 
nonvanishing only on the hypersurface in phase-space defined by V(x,p) = 0. 

In order to implement proper selfenergies for the nucleons in the p production 
processes we follow ref. [[24J and separate the mean-fields into a local part and an 
explicit momentum-dependent part, i.e. 



U s (x,p) = U«(x) + Ur(x,P) 



U,(x,p) = U»(x)+UM D (x,p). (6) 
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where the local mean-fields are determined by the usual Hartree equation: 

u si x ) = -9sVh(x) 

U%{x) = g v rf(x) (7) 

with 

m 2 s a H (x) + B s a 2 H (x) + C s a 3 H (x) = g s p s (x) 

™ 2 srf{x) =g v j tl (x). (8) 

In the above equations the scalar density p s (x) and the current are given in 

terms of the phase-space distribution function by 

4 r 



(27T) 



4 r 

3n{x) = J d 4 p U fl (x,p)f(x,p). (9) 

The free parameters of the above expressions are fixed to reproduce the sat- 
uration properties of nuclear matter, the empirical proton-nucleus optical potential 
as well as the density dependence of U s and from Dirac-Brueckner theory f26f . 



Explicit values for these parameters are given in p3| , p4| , g5j . The actual results for 



U s (p) and the zero'th component of the vector field Uq(p) are displayed in Fig.l for 
p (~ 0.17/m- 3 ),2p , and 3p . 

As mentioned before, the phase-space distribution function for the antiprotons 
fp(x,p) is assumed to follow an equation of motion as in (1), however, with scalar 
and vector potentials of different strength, i.e. 

Uf(x) = -g p s a H (x) 

K(x) = £uf(x) (10) 

where g& and g% are treated as independent parameters and are fixed in comparison 
to the experimental data (see below). The collision term (r.h.s. of eq. (1)) for the 
antiproton phase-space distribution besides elastic scattering also includes a direct 
coupling to the nucleons which describes the p-annihilation channel (cf. discussion 
below) . 

The antiproton invariant differential multiplicity is obtained by summing inco- 
herently the elementary antiproton multiplicities over all collisions and integrating 
over the residual degrees of freedom 0. If we consider that the antiprotons are 
produced via processes of the type BB — > pp + iViV = p + 3 + 4 + 5 (B stands for 



either nucleon or A) we can write the antiproton invariant multiplicity as [17 

d 3 N{b) 
P d 3 pp 

1 , d l2 <7BB->p{\fs) 



[l-f(r,p' 3 ;t)] [l-f(rJU;t)} [l ~ f(r,p' 5 ;t) 
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Here, the quantities p\ denote baryon momenta in the individual BB center-of-mass 
system which have to be transformed into the laboratory frame or the midrapidity 
frame, respectively, and s = {P1+P2) 2 is the squared invariant energy of the collision. 
Finally, the antiproton invariant differential cross section is obtained by integrating 
the differential multiplicity (|TTD over the impact parameter b. 

In order to proceed in the evaluation of QTTp we assume, as in refs. |Ll], [15], [17] ], 
that the differential elementary antiproton cross section is proportional to the phase- 
space available for the final state: 

E 3 E 4 E 5 E dl2(TBB ^P^ s ) 
p d 3 p p d 3 p 3 d 3 p 4: d 3 p 5 

= (TBB^p(Vs) - ~ — (5 4 (pi +P2-P3-P4-P5- Pp) ■ (12) 

16 R^Ws) 



The factor l/i? 4 (y / s) is the 4-body phase-space integral |3(J and has been 
included in order to normalize the differential distribution. It should be noted that 
the 4-body phase-space integral contrary to ref. [H] now strongly depends on the 



quasi-particle properties in the medium defined by eqs. (2) - (5) and (10). The 
total cross section for antiproton production UBB~*p{\fs) has been extracted from the 
experimental data corresponding to the inclusive process pp — > p+X. Unfortunately, 
there are no data available at y/s — 4m < 1 GeV so that we perform an extrapolation 



to lower energies. In line with [17j, [22], [Tj| we adopt the parametrization: 



a 



pp- 



,p+x(V~s) = 0.01 (yG - v ^) 1 - 846 [mb] (13) 

and assume that a pp ^p + x(yfs) = ubb^p p +nn(\^)- Whereas s/so = 3.7532 GeV for 
free particles, we replace ^/sq by the corresponding threshold ^Jso(p,pi) for the quasi- 
particles as defined by eqs. (2) - (5) in the medium. Though the latter assumption is 
not controlled by experimental data so far - and can hardly be measured - it is well 
in line with the concept of a phase-space dominated elementary production cross 
section. 

Apart from the perturbative production scheme discussed above, the antipro- 
tons produced in individual baryon-baryon collisions can be annihilated on their 
way out of the dense nuclear medium into the continuum. Since a proper treatment 
of p-annihilation is decisive for a comparison with experimental data, we perform 
this task nonperturbatively in the following way: For each baryon-baryon collision 
event i at space-time position we evaluate the differential production probability 
Pi{pj) for an antiproton with momentum pj on a grid in momentum space. Then 
each grid point in momentum space j (for each collision event i) is represented by 
an antiproton testparticle which is propagated in time under the influence of its 
selfenergies (10) according to the respective equations of motions for point particles 
which provides individual trajectories ry(t). In order to account for annihilation, 
the individual probabilities Pi{pj] t) are integrated in time according to the following 
rate equation 

d 4 r 

—Pi(pft) = -j^y J d z pv 12 a ann {p h p) f{r ij (t),p]t)P l {pft) (14) 
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where vu is the Lorentz-invariant relative velocity, f(r,p;t) the baryon phase-space 
distribution from eq. (1) and a ann the p annihilation cross section which is taken 



from the experimental data |31] and parametrized as a function of \fs — y/so- Though 
it will remain a matter of debate, how cr ann might change in the nuclear medium, we 
follow our concepts above and use the free annihilation cross section, however, mod- 
ify -^/so in the medium according to the quasi-particle properties described above. 
We note that the validity of this concept can in part be controlled via the exper- 
imental mass dependence of the p yield when comparing light and heavy systems. 



For further details we refer the reader to ref. p2 |. 

We have applied the above mentioned formalism to evaluate the antiproton 
cross section for the reactions p + 12 C and p + 63 Cu at bombarding energies of 
5, 4, and 3.5 GeV. The corresponding invariant cross sections in comparison with 
the data of ref. Jl0[ are shown in Fig. 2 (a + b) as a function of the momentum 



of the emitted antiproton in the lab frame at = 0°, assuming free antiprotons, 

1. e. gP, g% = 0. The calculations slightly underestimate the experimental data, but 
already approximately reproduce the shape of the momentum-spectra as well as the 
dependence on bombarding energy and mass. When adopting a slightly attractive 
scalar p selfenergy of - 50 to - 100 MeV the reproduction of the data improves at 
all energies significantly which is exemplified for 4.0 GeV by the dashed line in Fig. 

2. We note that in the above comparison we cannot distinguish between scalar and 
vector antiproton selfenergies because both yield similar results for the p spectrum 
if the same Schrodinger-equivalent optical potential is achieved. Furthermore, when 
using antiproton selfenergies in line with the relativistic mean-field theories i.e. 
changing only the sign of nucleon vector potential, we overestimate the p yield by 
more than an order of magnitude at all energies for both systems. 

We now turn to the nucleus-nucleus case. The calculated antiproton invariant 
differential cross section for the reaction 28 Si+ 28 Si at 2.1 GeV/A and Ni + Ni at 1.85 
GeV/u is shown in Figure 3 in comparison to the experimental data of ref.[|8| and 
ref. |ll| ; the upper lines represent the results of the calculations for free antiprotons 
without including any reabsorption. When taking care of antiproton annihilation 
according to eq. (14) the yields drop to the lower full lines which now underestimate 
the data sizeably. However, using attractive scalar (or vector) selfenergies at p = p 
of about - 100 to -150 MeV we nicely reproduce the data again. Since the two 
systems studied differ quite substantially in mass we infer that the description of p 
annihilation appears to be sufficiently accurate. 

We note in passing that the contribution due to collisions involving resonances 
is a factor 10 larger than those involving only nucleons thus confirming again the 



main results of ref. |L2|, [L7j . A detailed investigation of the baryonic decomposition 
for the present reactions is given in ref. |32] . 

The different value for the attractive antiproton field at p = po in p + A 
and A + A reactions is due to the fact that in p + A collisions the antiprotons 
move with momenta of 1 - 2 GeV/c with respect to the nuclear medium, whereas 
in A + A collisions the antiprotons are almost at rest in the nucleus-nucleus center- 
of-mass frame. In view of uncertainties of our present studies with respect to the 
elementary production cross sections close to the thresholds we provide areas for 
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the antiproton Schrodinger equivalent potential at p = po in Fig. 4, as extracted 
from the comparison with the experimental data for p + A and A + A reactions. 
These areas are far from the values expected by charge conjugation from the familiar 
a — co model [[Z7| (dashed line) and thus exclude relativistic mean-field models with 
the same parameter-sets for nucleons and antinucleons. However, our extrapolated 
values are well in line with a Schrodinger-equivalent potential (solid line in Fig. 4) 
as extracted from the dispersion relation (P: principle value) 

whereas the imaginary part is determined from the annihilation rate at density p 
according to 

2Im(U p (E,p )) = -| a ann (E) p (16) 
with E = y/p' 2 + m' 2 and a ann from ref. |31 1- The real part of the p-potential (as 



shown in Fig. 4 by the solid line) is well in line with the potential analysis for p + 
A reactions [3^ . 



In summary, we have evaluated the differential cross section for p production 
for p-nucleus and nucleus-nucleus reactions in the subthreshold regime by consider- 
ing incoherent on-shell baryon-baryon production processes involving nucleons and 
A's with their in-medium quasi-particle properties and treating p annhilation non- 
perturbatively. The quasi-particle properties of the nucleons are fixed within our 
study by the nuclear saturation properties, the proton-nucleus empirical potential as 
well as Dirac-Brueckner calculations at higher density. In comparing our calculations 
to the most recent data from KEK [10| and SIS [1 j] we find a consistent description 



of the experimental results with a rather weak attractive potential for the antipro- 
tons which is almost perfectly in line with a dispersive potential extracted from the 
dominant imaginary part of the antiproton selfenergy due to annihilation. Though 
we cannot exclude uncertainties due to uncertain elementary production cross sec- 
tions we can infer, that relativistic mean-field theories - which predict antiproton 
selfenergies according to charge conjugation (and violate unitary) - are definitely 
inadequate for the description of antimatter in a dense baryonic environment. 

We note in closing that the antiproton production studies at the AGS |34], 



35], |36] around 15 GeV/u might yield further information on the dynamics and 
selfenergies of antiprotons at even higher densities although these processes occur 
far above threshold. 



The authors acknowlegde valuable discussions with A. Gillitzer, P. Kienle, W. Konig 
and A. Schroter as well as their information on experimental results prior to publi- 
cation. 
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Figure Captions 

Fig.l: Scalar and vector selfenergies U s {p) and U (p) for nucleons at different 
densities in units of po ~ 0.17/m -3 as used in the relativistic transport 
equation (1). 

Fig. 2: Invariant cross section for antiproton production in the reactions p+ 12 C 
and p+ 63 Cu at = 0° as a function of the antiproton momentum p in 
the lab system. The experimental data are taken from ref. |Hj] and 
correspond to bombarding energies of 5.0 GeV , 4.0 GeV and 3.5 GeV. 
The full lines represent calculations for free antiprotons. The dashed 
lines indicate the result for an antiproton selfenergy of - 100 MeV at 4.0 
GeV. 

Fig. 3: Invariant cross section for antiproton production in the reaction 28 Si+ 28 Si 
at 2.1 GeV/u and Ni + Ni at 1.85 GeV/u for = 0° as a function of 
the momentum p of the emitted antiproton in the lab-system. The ex- 



perimental data have been taken from refs. || and |TT| , respectively. The 
upper lines indicate the calculated cross section for free antiprotons with- 
out reabsorption whereas the lower solid line is obtained when including 
p annihilation. The dashed line represents the cross section adopting an 
attractive potential of the antiproton of - 150 MeV. 

Fig. 4: Comparison of our extracted values for the Schrodinger equivalent an- 
tiproton potential from p + A and A + A reactions with the prediction 
from the a— u model (dashed line) and the dispersive potential according 
to eq. (15) (solid line). 
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